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SUMMARY 

'  This  Note  presents  an  attempt  to  develop  a  rotor  wake 
requires  only  modest  computing  time  and,  therefore,  is  sui.jt/r 
Free  and  rigid  wake  models  are  discussed,  together  with  vni  ,ou 
in  order  to  reduce  the  computation  time  to  acceptable  lev  Is  for  simulation  studies.  Using 
a  tip  vortex  wake  model  and  the  Biot-Savart  relation,  induced  velocities  are  calculated  in 
the  rotor  plane.  It  is  then  assumed  that  the  tip  vortex  wake  takes  the  form  of  a  skewed 
helix.  Such  a  simple  prescribed  model  requires  little  computing  time  and  the  accuracy  of 
its  results  when  compared  with  more  elaborate  free  wake  models  shows  considerable 
promise. 
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NOTATION 

Core  size  of  the  vortex  element  ij 
Length  of  the  vortex  element  ij 
Radial  distance  from  the  rotor  hub 

Distances  of  the  point  P  from  the  end  points  of  the  vortex  element  ij  as  shown  in 
Figure  6 

Outer  radius  of  the  rotor 
Induced  velocity 

Velocity  components  in  x,  y,  and  z  direction  respectively 
The  cartesian  co-ordinate  system 
Co-ordinates  of  the  point  P 

Tip-path-plane  angle 
Vortex  strength  of  the  element  ij 
Advance  ratio 
Azimuth  angle 

Rotational  speed  of  the  rotor  blade 
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1.  INTRODUCTION 

In  flight  simulation  studies,  the  forces  and  moments  acting  on  a  helicopter  rotor  blade  are 
often  estimated  on  the  assumption  that  the  distribution  of  induced  velocity  over  the  disc  is 
uniform.  Experimental  studies  have  revealed,  however,  that  the  assumption  is  not  well  founded, 
especially  in  forward  flight.  The  only  justification  for  using  the  uniform  inflow  distribution  is 
that  the  calculations  are  simple  and  the  results  obtained  using  this  assumption  are  in  fair  agree¬ 
ment  with  experimental  data.  Since  the  theory  gives  reasonable  average  rotor  thrust,  blade 
loads  and  moments,  it  can  be  used  satisfactorily  as  a  first  approximation  in  rotor  performance 
and  flight  simulation  calculations.  However,  to  obtain  a  more  fundamental  picture  of  the 
mechanism  of  rotor  induced  flow,  analysis  of  the  rotor  vortex  wake1  must  be  considered.  Though 
it  has  been  recognized  for  many  years  that  the  detailed  structure  of  the  rotor  downwash  could 
have  a  substantial  influence  on  rotor  characteristics  such  as  performance  and  blade  flapping 
motion,  early  efforts  at  computation  of  the  induced  flow  in  forward  flight  were  limited  in  utility 
because  it  was  impracticable  to  account  for  individual  blades.  The  advent  of  the  high-speed 
electronic  digital  computer  has  made  possible  the  straightforward  approach  of  tracing  the  vortex 
filaments  trailed  by  each  blade,  and  using  the  Biot-Savart  relation  to  obtain  velocities. 

In  the  following  sections  wake  models,  both  free  and  prescribed,  are  discussed,  together 
with  various  simplifying  assumptions  to  reduce  the  computational  time  required.  The  computa¬ 
tional  procedure  to  calculate  the  induced  velocity  due  to  vortex  elements  is  also  discussed.  It  is 
concluded  that  even  with  modern  fast  computers  the  computational  time  greatly  restricts  the 
application  of  the  free  wake  model.  On  the  other  hand  it  is  shown  that  a  simple  tip  vortex 
model  which  is  assumed  to  take  the  form  of  a  skewed  helix  gives  satisfactory  accuracy,  compared 
with  the  more  complex  free  wake  model.  It  is  considered  that  a  simple  prescribed  rotor  wake 
model  for  forward  flight  is  substantially  in  accordance  with  the  true  state  of  the  flow  generated 
by  a  rotor  blade  and  allows  downwash  calculations,  with  acceptable  computing  time,  for  flight 
simulation  studies. 

2.  THE  VORTEX  WAKE  MODEL  FOR  FORWARD  FLIGHT 

Consider  the  rotor  blade  in  steady  forward  flight.1  Because  the  circulation  generally  varies 
along  the  span,  vortex  lines  leave  the  trailing  edge  and  spiral  downwards  beneath  the  rotor.  In 
addition,  however,  the  blade  experiences  oscillatory  aerodynamic  effects  even  in  steady  forward 
flight.  This  is  fundamental  for  rotary  wings  since  variations  in  relative  airspeed  occur  at  each 
blade  section  as  the  blade  rotates.  Flapping  and  blade  pitch  angles  are  made  to  vary  with  azimuth 
in  order  to  tilt  the  rotor  thrust  vector  in  the  desired  direction.  Hence  the  blade  incidence  varies 
as  a  function  of  rotor  azimuth,  creating  shed  vortices  in  the  wake.  Thus  as  the  blades  rotate,  a 
continuous  sheet  of  shed  and  trailing  vorticity  streams  from  each  section  of  the  blade  as  shown 
in  Figure  1 .  Similar  wake  structures  exist  for  other  blades,  with  the  aggregate  forming  the  com¬ 
plete  wake.  This  system  of  vortex  wake  is  the  one  which  makes  the  helicopter  aerodynamic 
analysis  so  much  more  complicated  than  the  fixed  wing  aerodynamic  problem.  Whereas  for  the 
fixed  wing  the  wake  is  assumed  to  lie  in  the  same  plane  as  the  wing,  the  helicopter  wake  is  blown 
below  the  plane  of  the  rotor.  The  exact  configuration  of  the  vortex  wake  depends  on  the  local 
velocity,  which  is  the  sum  of  the  induced  velocity  and  the  main  flow  due  to  helicopter  motion. 
Since  the  resultant  local  velocity  is  determined  in  part  by  the  spatial  distribution  of  the  wake 
vorticity,  a  closed  form  analytic  solution  for  both  the  downwash  and  wake  vorticity  is  virtually 
impossible.  Hence  numerical  methods  must  be  used  to  solve  the  problem. 

A  detailed  wake  model  for  calculating  the  induced  velocity  would  represent  the  rotor  blades 
by  lifting  surfaces  and  the  rotor  wake  by  vortex  sheets.  The  calculation  of  wake  geometry  would 
involve  the  computation  of  the  distortion  and  roll-up  of  these  vortex  sheets  due  to  their  own 
induced  velocities  and  those  of  the  lifting  surfaces.  As  the  vortex  sheets  roll  up  into  line  vortices, 
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viscous  effects  would  become  important  in  the  vortex  core.2  Because  the  above  model  would 
require  an  excessive  amount  of  CPU  time,  a  simplified  model  must  be  constructed.  For  mathe¬ 
matical  analysis,  the  rotor  blades  are  represented  by  lifting  lines  (bound  vortices),  and  the 
vortex  wake  (shown  in  Fig.  1)  is  idealized  by  a  set  of  trailing  and  shed  line  vortices  as  shown  in 
Figure  2.  The  trailing  vortex  lines  are  trailed  aft  along  streamlines  relative  to  the  rotor  blade 
while  the  shed  vortex  lines  are  shed  parallel  to  the  instantaneous  position  of  the  rotor  blade. 
These  vortex  filaments  are  free  to  convect  at  the  local  velocity,  which  is  the  sum  of  the  free 
stream  velocity  and  the  velocity  induced  by  the  trailing,  shed,  and  bound  vortices.  This  form  of 
rotor  wake  representation  is  one  of  the  most  complex.  The  computer  requirements  and  cost  to 
run  such  a  wake  model  are  normally  prohibitive.  This  point  will  be  very  clear  later  when  some 
results  for  a  much  simplified  version  of  this  model  are  presented.  Comparisons  of  entire  wake 
flows  predicted  using  this  (refer  Fig.  2)  wake  model  with  experimental  measurements  have 
indicated  that  use  of  a  coarse  mesh  resulted  in  poor  induced  velocity  predictions,  and  that  use 
of  a  fine  mesh  increased  running  time  to  an  unacceptable  level.3'4  Therefore,  the  full  mesh  wake 
was  used4  to  represent  the  wake  immediately  behind  the  blades,  and  a  modified  wake  model 
was  used  in  the  remainder  of  the  wake,  as  shown  in  Figure  3.  The  modified  wake  consists  of 
trailing  vortices  only.  It  is  known  from  experiments  that  the  strength  of  the  shed  vorticity  is 
small  in  comparison  to  the  trailing  vorticity.  Hence  the  influence  of  the  shed  vorticity  on  the 
wake  distortions  is  neglected.  The  locations,  circulations,  and  core  sizes  of  the  trailing  vortices 
in  the  modified  wake  are  determined  in  terms  of  the  final  values  of  the  full  mesh  wake.  The  wake 
induced  velocities,  wake  distortions,  and  other  calculations  are  essentially  the  same  for  both  the 
full  mesh  wake  and  for  the  modified  wake  portions  of  the  wake  model.  Piziali5  used  the  above 
model  in  his  rigid  wake  analysis.  Clark  and  Leiper,8  Landgrebe,4  and  Sadler3  used  a  similar  vortex 
model  in  their  free  wake  analyses.  Although  this  wake  model,  when  used  in  free  wake  analysis, 
gives  results  which  show  good  agreement  with  experimental  results,  the  large  computational  time 
required  limits  its  use,  thus  necessitating  further  simplification. 

A  further  simplification  of  the  rotor  wake  geometry  is  valid  if  certain  approximations 
based  on  experimental  results  are  introduced.  It  is  known  that  the  vorticity  trailed  by  a  rotor 
blade  tends  to  be  concentrated  towards  the  blade  tip  and  rolls  up  into  a  core  of  rotating  air 
(see  Ref.  17).  Most  of  the  contribution  to  induced  velocity  and  blade  loading  comes  from  this 
tip  vortex.  Representation  of  the  real  wake  by  a  tip  vortex  alone  is  a  fairly  good  approximation.7 
This  wake  model  is  shown  in  Figure  4.  Crimi7  and  Scully*  used  similar  idealized  wake  models. 
In  Figure  4,  it  is  indicated  that  the  vorticity  which  rolls  up  into  the  tip  vortex  must  clearly  be 
shed  as  a  sheet;  however,  little  is  known  about  the  roll-up  process.  Flow  visualization  pictures 
suggest  that  roll-up  occurs  close  to  the  blade.  In  this  paper  it  is  assumed  that  the  roll-up  takes 
place  immediately  behind  the  blade  trailing  edge,  that  is,  it  is  assumed  that  the  vortex  sheet 
roll-up  angle  (see  Fig.  4)  is  zero. 

3.  COMPUTATIONAL  PROCEDURE  TO  CALCULATE  THE  INDUCED 
VELOCITY  DUE  TO  VORTEX  ELEMENTS 

In  the  previous  section,  a  simplified  model  was  evolved  for  the  vortex  wake  of  a  rotor  blade. 
The  next  task  is  to  determine  the  geometrical  distribution  of  the  vorticity  in  the  flow  field.  The 
wake  geometry  is  calculated  by  an  iterative  numerical  process.  As  a  part  of  the  computational 
procedure,  the  tip  vortex  behind  each  blade  is  broken  into  convenient  straight  line  segments 
as  shown  in  Figure  5.  These  segments  were  chosen  to  be  sufficiently  small  so  that  for  purposes 
of  computation  of  the  wake-induced  velocities  they  may  be  considered  as  rectilinear  vortices 
having  constant  circulation  along  their  length.  The  wake  configuration  at  any  instant  is  then 
defined  by  the  location  of  these  line  segments.  Each  vortex  filament  has  length  Uu  core  radius 
ai/,  and  vortex  strength  rw.  The  length  of  each  vortex  element  depends  on  how  finely  the  tip 
vortex  is  divided.  The  wake  and  bound  vortices  are  assumed  to  have  finite-sized  cores  of  rota¬ 
tional  fluid.  Various  sources7'*  estimate  the  vortex  core  radius  and  give  a  range  of  values. 
However,  for  the  case  of  a  rotor  in  forward  flight,  the  results  turn  out  to  be  very  insensitive  to 
vortex  core  size.  After  extensive  calculations,  Crimi7  assigns  a  representative  value  of  0  05  to 
the  ratio  of  core  radius  to  the  rotor  radius.  The  vortex  strength  F (f  of  the  filament  is  taken  as  the 
value  of  the  bound  vorticity  when  it  was  generated.  Once  formed,  the  strength  of  the  vortex 
elements  are  assumed  to  remain  unchanged  as  the  wake  develops.  Knowing  the  vortex  filament 


2 


length,  core  radius,  and  strength,  the  Biot-Savart  relation  can  be  used  to  calculate  the  induced 
velocity  at  any  point.  Figure  6  shows  a  typical  vortex  element  between  the  points  /  and  j.  Using 
cartesian  co-ordinates,  the  induced  velocity  components  due  to  this  vortex  filament  are  denoted 
by  v„  v„,  and  v2  in  the  x,  y  and  z  direction  respectively.  The  method  for  calculating  these  induced 
velocities  using  the  Biot-Savart  relation  is  well  known.3'4'7  For  completeness,  the  formulae  are 
given  in  Appendix  A.  However,  these  formulae  are  not  useful  when  it  is  necessary  to  calculate 
the  induced  velocity  at  certain  important  points  in  the  flow  field. 

For  the  case  of  calculating  the  velocity  induced  at  a  point  on  the  tip  vortex  by  the  two 
immediately  adjacent  vortex  segments,  the  use  of  the  straight  line  vortex  segments  is  inadequate, 
since  it  would  give  zero  induced  velocity  at  the  point  due  to  these  two  vortex  segments.  The  two 
vortex  segments  adjacent  to  the  point  are  therefore  replaced  by  a  circular  arc  vortex  segment 
passing  through  the  three  points  determined  by  the  ends  of  the  vortex  segments.  The  velocity 
induced  at  a  point  by  the  circular  arc  vortex  segment  is  given  in  References  3  and  7. 

Other  potentially  awkward  aspects  of  using  concentrated  line  vortices  arise  when  it  is 
attempted  to  compute  the  induced  velocity  very  close  to  a  vortex  filament.  For  this  case,  the 
effects  of  a  finite  vortex  core,  where  viscosity  becomes  important  and  the  Biot-Savart  relation 
no  longer  applies  directly  should  be  included.  A  useful  first  approximation  to  the  induced 
velocity  inside  a  vortex  core  is  to  assume  solid  body  rotation  of  the  vortex  core  with  an  induced 
velocity  at  the  outer  edge  of  the  vortex  core  equal  to  the  result  given  by  the  Biot-Savart  relation 
at  that  point.3  Thus  for  most  cases  it  isonly  necessary  to  check  whether  or  not  the  point  at  which 
the  induced  velocity  must  be  calculated  is  inside  the  vortex  core,  and  apply  either  the  Biot- 
Savart  relation  or  the  solid  body  rotation  model. 

4.  FREE  WAKE  METHOD 

In  the  free  wake  method  the  computations  are  initiated  by  specifying  the  initial  wake 
configuration.  The  velocity  contributions  of  all  vortex  elements  are  calculated  at  each  reference 
point  (end  points  of  vortex  filaments).  Then  these  points  are  allowed  to  propagate  with  the 
computed  velocity  generating  a  new  vortex  geometry.  This  process  is  repeated  until  a  wake  is 
generated  which  is  unchanged  from  one  iteration  to  the  next.  This  wake  is  taken  as  the  equi¬ 
librium  form  and  is  used  in  rotor  downwash  calculations. 

As  an  example  the  rotor  wake  for  a  four-bJaded  Wessex  rotor  is  shown  in  Figure  7.  The 
important  parameters  for  the  rotor  are: 

Advance  ratio  0-215 

Rotor  forward  tilt  3  •  5  deg 

Rotor  speed  22-2  rad/s 

Weight  53  400  N 

Rotor  radius  3  ■  5  m 

The  iterative  vortex  wake  solution  converges  for  this  case.  It  takes  36  minutes  CPU  time 
on  the  PDP  DEC  10  system.  For  clarity,  the  tip  vortex  of  only  one  of  the  blades  is  presented 
in  Figure  7.  The  wake  elements  become  distorted  as  a  result  of  interactions  with  other  elements 
in  the  wake.  The  large  displacements  of  the  wake  elements  two  to  three  rotor  radii  downstream 
must  be  viewed  with  caution  considering  stability  and  viscous  dissipation  problems  involved 
in  that  region  and  beyond.  These  large  downstream  distortions,  however,  do  not  significantly 
influence  the  flow  field  at  the  rotor  where  rotor  downwash  calculations  are  needed  in  perfor¬ 
mance  computations. 

Using  the  converged  wake,  the  induced  velocity  can  be  calculated  at  any  reference  point, 
again  applying  the  Biot-Savart  relation.  Time  averaged  induced  velocities  at  various  points  in 
the  rotor  plane  were  calculated  for  this  case  (refer  Fig.  11).  These  results  are  discussed  in  detail 
in  succeeding  sections. 

In  the  free  wake  method,  the  computational  procedure  is  fundamentally  an  iterative  process. 
Induced  velocities  at  wake  reference  points  are  calculated  using  the  previously  calculated  wake 
geometry  and  then  these  points  are  allowed  to  propagate  to  constitute  the  new  wake  geometry. 
The  process  is  repeated  until  convergence. 

The  converged  wake  is  the  one  which  exists  in  force-free  equilibrium*  and  the  aim  is  to  find 
the  arrangement  of  the  vortex  segments  in  space  which  will  satisfy  this  condition.  This  presents 
a  formidable  computational  problem  since  every  vor'cx  segment  is  influenced  by,  and  influences, 
every  other  segment. 
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Even  using  modem,  fast  computing  machines,  the  time  required  in  the  free  wake  analysis 
restricts  its  application.  For  example,  depending  on  the  advance  ratio,  4  to  10  iterations  may  be 
needed  in  the  rotor  wake  geometry  calculation.  This  means  30  to  50  minutes  CPU  time  on  the 
PDP  DEC  10  system.  The  lower  the  advance  ratio,  the  more  time  is  required  for  wake  vortex 
interaction.  Thus  in  rotor  downwash  calculations,  the  major  portion  of  the  computer  time  goes 
into  the  determination  of  the  wake  geometry. 


5.  RIGID  WAKE  METHOD 

If  the  wake  geometry  can  be  prescribed,  based  on  previous  theoretical  calculations  and 
experimental  result,  much  time  and  labour  can  be  saved.  This  is  particularly  necessary  where 
induced  downwash  calculations  form  a  part  of  the  overall  problem  of  helicopter  flight  simulation 
studies.  For  flight  simulation  the  overall  mathematical  model  is  divided  into  simple  realistic 
sections  (e.g.  rotor  aerodynamics,  fuselage  aerodynamics,  helicopter  control  system,  dynamics 
of  slung  bodies)  so  that  the  problem  can  be  accommodated  on  a  medium-sized  computer  and 
results  expected  in  a  reasonable  time.  With  this  application  in  view  the  various  simple  prescribed 
wake  methods  described  here  are  studied  and  the  accuracy  of  their  results  determined  by  com¬ 
paring  with  the  results  of  the  free  wake  model. 

The  concept  of  the  prescribed  wake  method  is  simple.  In  the  rotor  wake,  the  vortex  fila¬ 
ments  are  distributed  in  some  complex  form  as  shown  in  Figure  7.  This  has  to  be  arranged  into 
a  simple  logical  form  so  that  it  closely  approximates  the  actual  rotor  wake.  Willmer®  and  Loewy10 
used  infinite  line  vortices  to  represent  the  prescribed  wake.  Coleman  el  al.,n  and  Heyson  and 
Katzoff12  used  a  simple  pattern  consisting  of  circles  and  straight  lines.  Piziali,5  and  Kocurek 
and  Tangier13  used  a  method  in  which  the  vortex  filaments,  both  shed  and  trailing,  are  broken 
into  a  number  of  straight  elements,  so  that  the  wake  has  the  appearance  of  a  twisted  spider’s 
web.  Cook14  used  circular  vortex  rings  and  ring  segments  to  represent  the  vortex  wake  in  forward 
flight  downwash  calculations.  Recently,  Rayner15-18  used  circular  rings  to  study  the  hovering 
flight  of  birds  and  insects,  and  elliptical  rings  to  generate  the  wake  of  the  forward  flight  of 
animals.  This  interesting  work  could  be  classified  as  a  semi-rigid  wake  model. 

The  decomposition  of  the  spiral  trailing  vortex  filaments  into  simple  circular  elements  for 
the  helicopter  flight  may  be  easily  demonstrated.  In  Figure  8,  the  spiral  trailing  tip  vortex  fila¬ 
ments  are  traced  as  they  are  generated  by  a  four-bladed  advancing  rotor  (see  Ref.  18). 

Considering  the  first  section  of  tip  vortex  from  the  blade  2,  it  could  easily  be  approximated 
to  a  semi-circular  element.  Furthermore,  the  contributions  of  additional  spirals  may  be  approxi¬ 
mated  to  complete  circles.  On  this  basis.  Cook  developed  a  prescribed  model  as  shown  in 
Figure  9.  Results  obtained  by  Cook  using  this  model  are  also  presented  in  Figure  11.  Although 
the  geometric  differences  between  the  distorted  (Fig.  7)  and  undistorted  (Fig.  9)  wakes  are 
significant,  it  has  been  found1  that  the  differences  are  generally  unimportant  when  it  comes  to 
the  determination  of  rotor  downwash  for  conventional  forward  flight  conditions.  This  cannot 
be  generalized,  however,  to  include  flight  conditions  which  place  portions  of  the  wake  in  the 
plane  of  the  rotor,  for  example  during  manoeuvres. 

The  hovering  rotor  generates  approximately  a  quarter  to  a  third  of  its  lift  over  the  outer 
10%  of  the  blades.  The  passage  of  a  tip  vortex  close  to  the  blade  in  this  tip  region  has  been 
shown1  to  influence  significantly  the  blade  tip  loading.  In  forward  flight  with  forward  tilt  of  the 
rotor  disc,  the  blades  are  moving  away  from  the  wake,  and  therefore  are  less  influenced  by  it. 
The  vortices  do  not  pass  continuously  under  the  blades  in  the  predominant  loading  regions.  In 
fact,  for  rotors  with  a  small  number  of  blades  flying  at  high  advance  ratios,  there  are  no  blade- 
vortex  interactions  over  a  large  range  of  the  blade  azimuth  travel.  This  makes  the  rigid  wake 
model  a  fair  proposition  in  forward  flight  simulation. 

As  stated  earlier,  the  main  aim  of  the  work  reported  here  is  to  develop  a  rotor  wake  model 
so  that  downwash  calculations  can  be  made  in  a  few  minutes  of  computer  time.  If  interest  is 
limited  to  calculating  the  velocity  field  in  the  rotor  plane  a  simple  vortex  wake  geometry,  in  the 
form  of  a  skewed  helix,  can  be  used  in  combination  with  the  Biot-Savart  relation  to  produce  a 
satisfactory,  simple  model.  The  results  from  this  model  are  compared  with  those  of  Cook  and 
free  wake  calculations  to  assess  their  accuracy.  From  Figures  7  and  10  it  is  clear  that  the  skewed 
helix  closely  represents  the  free  wake  near  the  rotor  but  that  it  differs  increasingly,  from  the 


free  wake,  with  distance  below  the  rotor  plane.  However,  this  divergence  is  not  important  since 
the  vortex  elements  far  away  from  the  rotor  contribute  very  little  to  the  induced  velocity.  The 
reason  for  this  is  that  the  induced  velocity  calculated  at  a  point  by  the  Biot-Savart  relation  is 
inversely  proportional  to  its  distance  from  the  vortex  element.  The  formulae  for  calculating 
the  wake  geometry  are  given  in  Appendix  B.  Induced  velocities  calculated  using  the  prescribed 
wake  model  are  also  presented  in  Figure  11,  for  the  Wessex  rotor,  again  for  the  same  operating 
conditions. 

Figure  1 1  shows  the  time-averaged  induced  velocities  in  the  rotor  plane  calculated  using 
free  vortex  and  prescribed  helical  wake  models,  along  with  Cook's  circular  ring  vortex  model. 
Glauert's  classical  trapezoidal  downwash  distribution  is  also  presented  in  the  figur».  Induced 
velocity  as  a  function  of  blade  span  is  shown  at  12  azimuth  stations.  The  graphs  illustrate  the 
magnitude  of  the  difference  between  various  approaches  and  show  that  the  basic  agreement 
between  the  free  vortex  and  prescribed  helical  wake  induced  velocity  results  is  fair.  The  charac¬ 
teristic  upwash  at  the  front  of  the  rotor  disc  is  well  defined.  By  far  the  largest  difference  between 
the  free  vortex  and  prescribed  helical  wake  results  is  at  an  azimuth  angle  of  300°. 

With  Cook’s  vortex  ring  model,  over  the  front  half  of  the  rotor  disc,  the  amplitude  of  the 
induced  velocity  fluctuates  along  the  blade  span.  These  fluctuations  are,  however,  not  observed 
in  the  results  calculated  here  using  the  free  vortex  and  prescribed  helical  wake  models. 

6.  CONCLUDING  REMARKS 

A  number  of  mathematical  models  for  calculation  of  rotor  downwash  in  forward  flight  has 
been  investigated.  It  is  found  that  even  the  tip  vortex  model  in  free  wake  analysis  requires  about 
50  minutes  CPU  time  on  a  PDP  DEC  10  for  the  Wessex  helicopter.  For  flight  simulation  studies 
this  amount  of  computer  time  is  unacceptable.  Various  ways  to  reduce  the  time  are  being 
investigated.  The  most  useful  approach  for  practical  applications  seems  to  be  an  empirical  one 
based  on  a  prescribed  wake  type  of  approximation.  The  development  of  a  satisfactory  method 
requires  the  knowledge  of  what  kind  of  wake  geometries  occur  under  various  flight  conditions. 
This  in  turn  can  only  be  acquired  by  the  extensive  analysis  of  rotor  wake  experimental  data, 
theoretical  results,  and  physical  reasoning.  Preliminary  results  carried  out  using  a  simple  skewed 
helix  model  are  very  encouraging. 
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APPENDIX  A 

Let  <  and  j  be  the  end  points  of  a  straight  line  vortex  filament  and  P  be  the  point  at  which 
the  induced  velocity  is  to  be  computed.  These  three  points  can  be  arbitrarily  oriented  in  space 
as  shown  in  Figure  6.  The  induced  velocity  components  vz.  »>.  and  rz  at  P  due  to  the  vortex 
element  I\/  are  given  by  the  Biot-Savart  relations:4,5 

p 

»’*  =  ~  c [(>'„  -  yt)(2p  -  z,)  -  (Zp  -  zi)(yp  -  j’/)] 

47T 

p 

»V  =  ~  C  [(zp  -  Zi)(xp  -  x ,)  -  ( xp  -  Xi)(zp  -  z/)] 

4  n 

p 

»•*  =  C[(xv  -  x,)(yp  -  yt)  -  (yp  -  yi)(xp  -  jtj)] 

47T 

where 

c=i _ ('<  +  n) _ 

R  rtrfcw  +  Ci  +  Cz  +  C3) 

Ci  =  (xp  -  x,)(xp  -  xj) 

Cz  =  (yp  -  yi)(yp  -  yi) 

C3  =  ( zp  -  2()(zp  -  2,) 

n  =  [(xp  -  *<)2  +  (yv  —  yt )2  +  ( zp  —  zt)2]1  2 

n  =  [(xp  -  x i)2  +  (yp  -  yi)2  +  (zp  -  z/)2]1'2 


APPENDIX  B 


Co-ordinates  of  the  skew  helix  wake  reference  point  Pt$  (refer  Fig.  5)  are  calculated  using 
the  following  formulae  :3,5,7 

x(i,j)  =  cos  jV  +  O'  -  1)  —  —  (»'  —  1)A^~|  -(-  (i  -  l)Ai/</i  cos  «t 
A‘J)  =  sin  jjfr  +  U  ~  1) 

{/Nb  w  v/21 

^in,T+^TiT2— —  j  J 

where  <l>  =  Azimuth  angle  as  shown  in  Figure  5 
Ab  =  Number  of  blades 
p.  =  Advance  ratio 

=  Blade  tip-path-plane  angle  as  shown  in  Figure  5 

Aip  =  Angle  (2 tr/N^);  Na  being  the  number  of  azimuth  stations  into  which  rotor  plane  is 
divided 

W  =  Weight  of  the  rotorcraft 
p  —  Air  density 
ft  =  Rotor  angular  speed 
R  =  Rotor  radius 
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FIG.  2  ROTOR  WAKE  STRUCTURE  -  IDEALIZED  LINE  VORTEX  MODEL 


FIG.  3  COMPLE 
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FIG.  6  FLOW  INDUCED  BY  A  VORTEX  ELEMENT 


FIG  7  DISTORTED  WAKE  BELOW  THE  ROTOR  AFTER  THE  CONVERGENCE 
OF  THE  FREE  WAKE  SOLUTION 
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